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The Monte Carlo simulation of N point vortices with square periodic boundary conditions is performed where N is order of 
100. The clustering property is examined by computing the L function familiar in the field of spatial ecology. The case of a 
positive value of L corresponds to the state of clustering and the Onsager's negative temperature. 
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1 Monte Carlo simulation of point vortices in a periodic box 
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The motion of point vortices (PVs) in a periodic box has been studied numerically in [1,2]. The K and L functions, which are 
familiar in the field of spatial ecology, are introduced in order to quantitate clustering of PVs. Periodic boundary conditions 
y—{ , guarantee spatial homogeneity and we need not to include the edge corrections in the K function. 

C*") ' In this report, the Monte Carlo simulation is performed for 100 identical PVs in order to obtain the probability distribution 
function (PDF) for the hamiltonian. Yatsuyanagi et al. [3] showed the PDF for two types of positive and negative PVs and the 
£^ ■ clustered state of the negative temperature. Figure 1 shows the histogram of the hamiltonian for TV identical PVs 
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c/3 ■ where Zi is the complex position of the PV, Re denotes the real part, a is the Weierstrass sigma function parametrized by two 

. ■ numrj ers wi = 1/2 and 0J2 = i/2, Q = w, N = 100, and the ensemble number is 10 4 . The il term denotes a rigid rotation 

C/3 ■ centered at each PV. Similarly to [3], the distribution with decaying tails is obtained. The figures 2-4 and 5-7 show the spatial 
distribution of PVs and the L functions for the minimum, median and maximum of the hamiltonian, respectively. The L and 

'" ^ ■ K functions are defined by 



L(r) = VK(r)/7r - r, (2) 
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K(r) = (\N)- l J2 E 0(r-\*-*i\), (3) 

. i=1 .7=1, 7^ 

, where A = N/ S, S is the area, and 9 is the step function. For complete spatial randomness (CSR), we have K = irr 2 and 
7-H ' L = 0. The positive (negative) values of L imply clustering (uniform spacing). 

, The right-side tail in Figure 1 corresponds to clustering of PVs, the positive value of L, and the state of the Onsager's 
' negative temperature. The left-side tail in Figure 1 shows uniform spacing of PVs, the negative value of L, and the positive 
J> , temperature. Figure 6 shows an oscillatory behavior of the L function for the median of H, where H = 2.61 is close to the 
' peak of the histogram (H = 2.613). Since there is a minimum in the hamiltonian for each pair of (i, j), the total hamiltonian 
. H is also lower bounded. However, we have no upper bound for H since a(z) ~ z for z ~ 0. In general, it is difficult to 
' judge whether the point distribution is clustered, completely spatially random, or uniformly spaced. The L function is found 
to be useful for the present purpose. 

The numerical simulation of 100 PVs with various conditions has been done in [1,2]. One of the typical behaviors observed 
in the case of positive and negative PVs with the same strength is the scattering and recoupling of the linearly moving pair 
of PVs. A gradient method for detecting stable stationary configurations of PVs, which corresponds to the system of point 
sources with a uniform sink, has been found by the author and will be reported elsewhere. 
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Fig. 1 Histogram of the hamiltonian of periodic point vortices with the Monte Carlo method. The right(left)-side tail denotes the negative 
(positive) temperature state. 
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Fig. 2 The distribution of PVs for the min- Fig. 3 The distribution of PVs for the me- Fig. 4 The distribution of PVs for the max- 
imum of the hamiltonian. dian of the hamiltonian. imum of the hamiltonian. 
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Fig. 5 The L function for the minimum of Fig. 6 The L function for the median of Fig. 7 The L function for the maximum of 
the hamiltonian. the hamiltonian. the hamiltonian. 
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